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Abstract 

The ANTARES telescope is well-suited to detect neutrinos produced in astrophys- 
ical transient sources as it can observe a full hemisphere of the sky at all times 
with a high duty cycle. Radio-loud active galactic nuclei with jets pointing almost 
directly towards the observer, the so-called blazars, are particularly attractive po- 
tential neutrino point sources. The all-sky monitor LAT on board the Fermi satellite 
probes the variability of any given gamma-ray bright blazar in the sky on time scales 
of hours to months. Assuming hadronic models, a strong correlation between the 
gamma-ray and the neutrino fluxes is expected. Selecting a narrow time window on 
the assumed neutrino production period can significantly reduce the background. 

An unbinned method based on the minimization of a likelihood ratio was applied 
to a subsample of data collected in 2008 (61 days live time). By searching for 
neutrinos during the high state periods of the AGN light curve, the sensitivity to 
these sources was improved by about a factor of two with respect to a standard time- 
integrated point source search. First results on the search for neutrinos associated 
with ten bright and variable Fermi sources are presented. 



Key words: ANTARES, Neutrino astronomy, Fermi LAT transient sources, 
time-dependent search, blazars 
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i 1 Introduction 



2 Neutrinos are unique messengers to study the high-energy universe as they 

3 are neutral and stable, interact weakly and therefore travel directly from their 

4 point of creation to the Earth without absorption. Neutrinos could play an 

5 important role in understanding the mechanisms of cosmic ray acceleration 
e and their detection from a cosmic source would be a direct evidence of the 
7 presence of hadronic acceleration. The production of high-energy neutrinos 
s has been proposed for several kinds of astrophysical sources, such as active 

9 galactic nuclei (AGN), gamma-ray bursters (GRB), supernova remnants and 

10 microquasars, in which the acceleration of hadrons may occur (see Ref. [1] for 

11 a review). 



12 Flat-Spectrum Radio Quasars (FSRQs) and BL Lacs, classified as AGN blazars, 

13 exhibit relativistic jets pointing almost directly towards the Earth and are 

14 some of the most violent variable high energy phenomena in the Universe [2] . 

15 These sources are among the most likely sources of the observed ultra high 

16 energy cosmic rays. Blazars typically display spectra with enhanced emission 

17 over two energy ranges: the IR/X-ray and MeV/TeV peaks. The lower energy 
is peak is generally agreed to be the product of synchrotron radiation from ac- 

19 celerated electrons. However, the origin of the higher energy peak remains to 

20 be clarified. In leptonic models [3], inverse Compton scattering of synchrotron 

21 photons (or other ambient photons) by accelerated electrons generates this 

22 high energy emission. In hadronic models [4] , MeV-TeV gamma-rays and high 

23 energy neutrinos are produced through hadronic interactions of the high en- 

24 ergy cosmic rays with radiation or gas clouds surrounding the source. In the 

25 latter scenario, a strong correlation between the gamma-ray and the neutrino 

26 fluxes is expected. The gamma-ray light curves of bright blazars measured by 

27 the LAT instrument on board the Fermi satellite reveal important time vari- 

28 ability on timescales of hours to several weeks, with intensities much larger 

29 than the typical flux of the source in its quiescent state [5]. 



30 This paper presents the results of the first time-dependent search for cosmic 

31 neutrino sources by the ANTARES telescope. The data sample used in this 

32 analysis and the comparison to Monte Carlo simulations are described in Sec- 

33 tion 2, together with a discussion on the systematic uncertainties. The point 

34 source search algorithm used in this time-dependent analysis is explained in 

35 Section 3. The search results are presented in Section 4 for ten selected can- 

36 didate sources. 
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2 ANTARES 



ss The ANTARES Collaboration completed the construction of a neutrino tele- 

39 scope in the Mediterranean Sea with the connection of its twelfth detector line 

40 in May 2008 [6]. The telescope is located 40 km off the Southern coast of France 

41 (42°48'N, 6°10'E) at a depth of 2475 m. It comprises a three-dimensional array 

42 of photomultipliers housed in glass spheres (optical modules [7]), distributed 

43 along twelve slender lines anchored at the sea bottom and kept taut by a buoy 

44 at the top. Each line is composed of 25 storeys of triplets of optical mod- 

45 ules (OMs), each housing one 10-inch photomultiplier. The lines are subject 

46 to the sea currents and can change shape and orientation. A positioning sys- 

47 tern based on hydrophones, compasses and tiltmeters is used to monitor the 

48 detector geometry with an accuracy of 10 cm. 

49 The main goal of the experiment is to search for high energy neutrinos with 

50 energies greater than 100 GeV by detecting muons produced by the neutrino 

51 charged current interaction in the vicinity of the detector. Due to the large 

52 background from downgoing atmospheric muons, the telescope is optimized 

53 for the detection of upgoing muons as only they can originate from neutrinos. 

54 Muons induce the emission of Cherenkov light in the sea water. The arrival 

55 time and intensity of the Cherenkov light on the OMs are digitized into hits 

56 and transmitted to shore. Events containing muons are selected from the con- 

57 tinuous deep sea optical backgrounds due to natural radioactivity and biolu- 

58 minescence. A detailed description of the detector and the data acquisition is 

59 given in [6,8]. 

eo The arrival times of the hits are calibrated as described in [9]. A LI hit is 

ei defined either as a high-charge hit, or as hits separated by less than 20 ns 

62 on OMs of the same storey. At least five LI hits are required throughout the 

63 detector within a time window of 2.2 /xs, with the relative photon arrival times 

64 being compatible with the light coming from a relativistic particle. Indepen- 

65 dently, events which have LI hits on two sets of adjacent or next-to- adjacent 
ee floors are also selected. 

67 The data used in this analysis were taken in the period from September 6 

es to December 31, 2008 (54720 to 54831 modified Julian days, MJD) with the 

69 twelve line detector. This period overlaps with the availability of the first 

70 data from the LAT instrument onboard the Fermi satellite. The corresponding 

71 effective live time is 60.8 days. Atmospheric neutrinos are the main source of 

72 background in the search for astrophysical neutrinos. These upgoing neutrinos 

73 are produced by the interaction of cosmic rays in the Earth's atmosphere. To 

74 account for this background, neutrino events were simulated according to the 

75 parametrization of the atmospheric neutrino flux from Ref. [10]. Only charged 
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76 current interactions of muon neutrinos and antineutrinos were considered. An 

77 additional source of background is due to downgoing atmospheric muons mis- 

78 reconstructed as upgoing. Downgoing atmospheric muons were simulated with 

79 the MUPAGE package [11]. In both cases, the Cherenkov light was propagated 
so taking into account light absorption and scattering in sea water [12]. 

si From the timing and position information of the hits, muon tracks are recon- 

82 structed using a multi-stage fitting procedure, based on Ref. [13]. The initial 

as fitting stages provide the hit selection and starting point for the final fit. The 

84 final stage consists of a maximum likelihood fit of the observed hit times and 

85 includes the contribution of optical background hits. 

se Upgoing tracks are also required to have a good reconstruction quality. The 

87 latter is quantified by a parameter, A which is based on the value of the 

as likelihood function obtained for the fitted muon (see Ref. [13] for details). The 

89 cumulative distribution of A for muons reconstructed as upgoing is shown 

90 in Figure 1 along with the simulated contributions from atmospheric muons 

91 and neutrinos. The angular uncertainty obtained from the muon track fit is 

92 required to be smaller than 1 degree. For this analysis, events are selected 

93 with A > —5.4. This value results in an optimal compromise between the 

94 atmospheric neutrino and muon background reduction and the efficiency of the 

95 cosmic neutrino signal with an assumed spectrum proportional to E~ 2 , where 

96 E u is the neutrino energy, which gives the best 5a discovery potential. The 

97 resulting sample consists of 628 events obtained in 60.8 days. The simulations 

98 indicate that the selected sample contains 60 % atmospheric neutrinos; the 

99 rest being mis-reconstructed atmospheric muons. 

100 The angular resolution of the reconstructed neutrino direction can not be 

101 determined directly from the data and has to be estimated from simulation. 

102 However, comparison of data and Monte Carlo in which the time accuracy of 

103 the hits was degraded by up to 3 ns constrains the uncertainty of the angular 

104 resolution to about 0.1° [14]. Figure 2 shows the cumulative distribution of the 

105 angular difference between the reconstructed muon direction and the neutrino 
we direction for an assumed spectrum proportional to E~ 2 . For the considered 

107 period, the median resolution is estimated to be 0.5 ± 0.1 degrees. 

108 The effective area for muon neutrinos is defined as the ratio between the rate 

109 of selected neutrino events and the cosmic neutrino flux. Figure 3 shows the 
no muon neutrino and antineutrino effective area of the ANTARES telescope as 
in a function of the declination of the source, after integrating over the energy 
n2 with an assumed spectrum proportional to E~ 2 between 10 GeV and 10 PeV. 
n3 In the flux limits (see Section 4), a conservative uncertainty on the detection 
H4 efficiency of about 30 % was taken into account. This number includes con- 
ns tributions on the uncertainty of the sea water optical parameters [12] and the 
lie OM properties such as efficiency and angular acceptance. 
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Track fit Quality (A) 



Fig. 1. Track fit quality (A) distribution for upgoing events in data (dots) and Monte 
Carlo samples (atmospheric muons: dashed line; atmospheric neutrinos: continuous 
line). Events are selected with an error estimate lower than 1 degree. The green 
dashed vertical line corresponds to the optimized event selection (A > —5.4). 
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Fig. 2. Cumulative distribution of the angle between the true Monte Carlo neutrino 
direction (a UMO ) and the reconstructed muon direction (a rec ) for an E~ 2 flux of 
upgoing neutrino events selected for this analysis. 

117 3 Time-Dependent Search Algorithm 

us The time- dependent point source analysis is performed using an unbinned 
n9 method based on a likelihood ratio maximization. The data are parametrized 
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Fig. 3. ANTARES muon neutrino and antineutrino effective area (continuous line) 
as a function of the declination of the source computed from the Monte Carlo 
simulation for an E~ 2 flux of upgoing muons selected for this analysis. The product 
of the effective area by the visibility (i.e. fraction of the time the source is visible 
at the ANTARES location) is shown with the dashed line. 

as a mixture of signal and background. The goal is to determine, at a given 
point in the sky and at a given time, the relative contribution of each com- 
ponent and to calculate the probability to have a signal above background 
in a given model. The likelihood ratio, A, is the logarithm of the ratio of the 
probability density for the hypothesis of signal and background (H sig+bkg ) over 
the probability density of only background (H b kg)'- 



N 



U ° 9 P{x t \H bkg ) 



jy Psigi^^iiti) (1 jy ~) Pbkg^ijti) 



bkg 



a 



126 where n S i g is the unknown number of signal events determined by the fit and 

127 N is the total number of events in the considered data sample. P s ig((Xi,ti) 

128 and Pbk g (5i,U) are the probability density functions (PDF) for signal and 

129 background respectively. For a given event i, U, Si and at represent the time 

130 of the event, its declination and the angular separation from the source under 

131 consideration. 



132 The probability densities P sig and P b k g are factorized into a purely directional 

133 and a purely time-related component. The shape of the time PDF for the 

134 signal event is extracted directly from the gamma-ray light curve assuming 

135 proportionality between the gamma-ray and the neutrino fluxes. It is assumed 

136 that the muon neutrino velocity in vacuum is equal to that of light in vacuum. 

137 For signal events, the directional PDF is described by the one dimensional 
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138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 



point spread function (PSF), which is the probability density of reconstruct- 
ing an event at an angular distance a from the true source position. The 
directional and time PDF for the background are derived from the data using 
the observed declination distribution of the selected events and the observed 
one-day binned time distribution of all the reconstructed muons respectively. 
Figure 4 shows the time distribution of all the reconstructed events and the 
selected upgoing events for this analysis. Once normalized to an integral equal 
to 1, the distribution for all reconstructed events is used directly as the time 
PDF for the background. Empty bins in the histograms correspond to periods 
with no data taking (i.e. detector in maintenance) or with very poor quality 
data (high bioluminescence or bad calibration). 
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Fig. 4. Time distribution of the reconstructed events. Upper histogram (black line): 
distribution of all reconstructed events. Bottom filled histogram (red): distribution 
of selected upgoing events. 

149 The statistical interpretation of the search result relies on simulated pseudo 

150 experiments (PE) in which the background events are randomly generated by 

151 sampling the declination and the time from the parametrization Pbkg{$i,ti) 

152 and the right ascension from a uniform distribution. Events from a neutrino 

153 point source are simulated by adding events around the desired coordinates 

154 according to the point spread function and the time distribution of the studied 

155 source. Systematic uncertainties (cf Section 2) are incorporated directly into 

156 the pseudo experiment generation. 

157 The null hypothesis corresponds to n sig = 0. The obtained value of Xdata on 

158 the data is then compared to the distribution of \{n sig = 0). Large values of 

159 ^data compared to the distribution of \{n sig = 0) reject the null hypothesis 

160 with a confidence level (C.L.) equal to the fraction of the number of PE above 
I" ^data- The fraction of PE for which \{n sig = 0) is above Xdata is referred to as 
162 the p-value. The discovery potential is then defined as the average number of 
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163 signal events required to achieve a p-value lower than 5a in 50 % of the PEs. 

164 In the same way, the sensitivity is defined as the average signal required to 

165 obtain a p-value less than that of the median of the X(n sig = 0) distribution 
lee in 90 % of the PEs. In the absence of evidence of a signal, an upper limit on 
167 the neutrino fluence is obtained and defined as the integral in energy and time 
lea of the flux upper limit with an assumed energy spectrum proportional to E~ 2 

169 from 10 GeV to 10 PeV. The limits are calculated according to the classical 

170 (frequentist) method for upper limits [15]. 

171 The performance of the time-dependent analysis was computed by applying 

172 this unbinned algorithm for a single source assuming a single square-shape flare 

173 with a width varying from 0.01 days to 84 days. The solid line in Figure 5 shows 

174 the average number of events required for a discovery from one source located 

175 at a declination of -40° as a function of the width of the flare. The numbers 

176 in the black line are compared to that obtained without using the timing 

177 information (dashed line). The flare timing information yields an improvement 

178 of the discovery potential by about a factor 2-3 with respect to a standard 

179 time-integrated point source search [14]. 
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Fig. 5. Average number of events (solid line) required for a 5a discovery (50 % 
probability) from a single source located at a declination of -40° as a function of the 
width of the flare period (at) for the 60.8 day analysis. These numbers are compared 
to that obtained without using the timing information (dashed line). 



wo 4 Search for Neutrino Emission from Gamma-Ray Flares 

181 The time-dependent analysis was applied to bright and variable Fermi blazar 

182 sources reported in the first-year Fermi LAT catalogue [16] and in the LBAS 

183 catalogue (LAT Bright AGN sample [17]). Sources were selected in the sky vis- 
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184 ible to ANTARES and that had at least one day binned gamma-ray flux in the 

185 high state periods greater than 80xl0~ 8 photons cm" 2 s" 1 above 100 MeV and 

186 showed significant time variability on time scales of days to weeks in the stud- 
is/ ied time period. A source is assumed variable in the LBAS catalogue when the 
las observation has a probability of less than 1 % of being a steady source. This list 
189 includes six flat-spectrum radio quasars and four BL-Lacs. Only four bright 
wo and nearby sources in the considered sample, PKS2155-304 [18], PKS1510- 

191 089 [19], 3C279 [20] and WComae [21], have been detected by the ground 

192 Cherenkov telescopes HESS, MAGIC or VERITAS. Table 1 lists the charac- 

193 teristics of the ten selected sources. 



AT 

JName 


UIGL name 


Class 


RA |°J 


Dec [ J 


Kedsmtt 


PKS0208-512 


J0210.8-5100 


FSRQ 


32.70 


-51.2 


1.003 


AO0235+164 


J0238.6+1636 


BLLac 


39.65 


16.61 


0.940 


PKS0454-234 


J0457.1-2325 


FSRQ 


74.28 


-23.43 


1.003 


OJ287 


J0855.4+2009 


BLLac 


133.85 


20.09 


0.306 


WComae 


J1221. 7+28.14 


BLLAc 


185.43 


28.14 


0.102 


3C273 


J1229. 1+0202 


FSRQ 


187.28 


2.05 


0.158 


3C279 


J1256.1-0548 


FSRQ 


194.03 


-5.8 


0.536 


PKS1510-089 


J1512.7-0905 


FSRQ 


228.18 


-9.09 


0.36 


3C454.3 


J2254.0+1609 


FSRQ 


343.50 


16.15 


0.859 


PKS2155-304 


J2158.8-3014 


BLLac 


329.70 


-30.24 


0.116 



Table 1 

List of bright variable Fermi blazars selected for this analysis [17]. 

194 The light curves published on the Fermi web page for the monitored sources [22] 

195 are used for this analysis. They correspond to the one-day binned time evolu- 

196 tion of the average gamma-ray flux above a threshold of 100 MeV since August 

197 2008. The high state periods are defined using a simple and robust method 

198 based on three main steps. Firstly, the baseline is determined with an iterative 

199 linear fit. After each fit, bins more than two sigma (o~bl) above the baseline 

200 (BL) are removed. Secondly, seeds for the high state periods are identified by 

201 searching for bins significantly above the baseline according to the criteria: 

(F — &f) > (BL + 2 * (?bl) and F > (BL + 3 * o BL ) (2) 

202 where F and op represent the flux and the uncertainty on this flux for each 

203 bin, respectively. For each seed, the adjacent bins for which the emission is 

204 compatible with the flare are added if they satisfy: (F — a F ) > (BL + obl)- 

205 Finally, an additional delay of 0.5 days is added before and after the flare in 
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206 order to take into account that the precise time of the flare is not known (1-day 

207 binned light curve). With this definition, a flare has a width of at least two 

208 days. Figure 6 shows the time distribution of the Fermi LAT gamma-ray light 

209 curve of 3C454.3 for almost two years of data and the corresponding selected 

210 high state periods. With the hypothesis that the neutrino emission follows the 

211 gamma-ray emission, the signal time PDF is simply the normalized light curve 

212 of only the high state periods. The third column of Table 2 lists the flaring 

213 periods for the ten sources found from September to December 2008. 




Fig. 6. Gamma-ray light curve (black points) of the blazar 3C454.3 measured by the 
LAT instrument onboard the Fermi satellite above 100 MeV for almost two years 
of data. The shaded histogram (blue) indicates the high state periods. The dashed 
line (red) represents the fitted baseline. 



214 The results of the search for coincidences between flares and neutrinos are 

215 listed in Table 2. For nine sources, no coincidences are found. For 3C279, a 

216 single high-energy neutrino event is found in coincidence during a large flare 

217 in November 2008. Figure 7 shows the time distribution of the Fermi gamma- 

218 ray light curve of 3C279 and the time of the coincident neutrino event. This 

219 event was reconstructed with 89 hits distributed on ten lines with a track 

220 fit quality A = —4.4. The particle track direction is reconstructed at 0.56° 

221 from the source location. The pre-trial p- value is 1.0 %. However, the post- 
222 trial probability computed taking into account the ten searches is 10 %; this 

223 occurrence is thus compatible with a background fluctuation. In the absence 

224 of a discovery, upper limits on the neutrino fluence were computed and are 

225 shown in the last column of Table 2. 
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Fig. 7. Gamma-ray light curve (dots) of the blazar 3C279 measured by the LAT 
instrument onboard the Fermi satellite above 100 MeV. The light shaded histogram 
(blue) indicates the high state periods. The dashed line (green) corresponds to the 
fitted baseline. The red histogram displays the time of the associated ANTARES 
neutrino event. 



226 5 Summary 



227 This paper presents the first time-dependent search for cosmic neutrinos using 

228 the data taken with the full twelve line ANTARES detector during the last four 

229 months of 2008. For variable sources, time-dependent point searches are much 

230 more sensitive than time-integrated searches due to the large reduction of the 

231 background. This search was applied to ten very bright and variable Fermi 

232 LAT blazars. One neutrino event was detected in time/direction coincidence 

233 with the gamma-ray emission in only one case, for a flare of 3C279 in November 

234 2008, with a post-trial probability of 10 %. Upper limits were obtained on the 

235 neutrino fluence for the ten selected sources. 
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Table 2 



Results of the search for neutrino emission in the ten selected sources. The meaning 
of the columns is the following: Vis: fraction of the time the source is visible at 
the ANTARES location; timePDF: high state periods of the light curve; LT: corre- 
sponding ANTARES live time in days; N(5<r): averaged number of events required 
for a 5(7 discovery (50 % probability); N^ a : number of observed events in time/angle 
coincidence with the gamma-ray emission. Fluence U.L.: Upper limit (90 % C.L.) 
on the neutrino fluence in GeV cm~ 2 . 



243 fur Bildung und Forschung (BMBF), Germany; Istituto Nazionale di Fisica 

244 Nucleare (INFN), Italy; Stichting voor Fundamenteel Onderzoek der Materie 

245 (FOM), Nederlandse organisatie voor Wetenschappelijk Onderzoek (NWO), 

246 the Netherlands; Council of the President of the Russian Federation for young 

247 scientists and leading scientific schools supporting grants, Russia; National 



14 



248 Authority for Scientific Research (ANCS), Romania; Ministerio de Ciencia 

249 e Innovation (MICINN), Prometeo of Generalitat Valenciana and MultiDark, 

250 Spain. We also acknowledge the technical support of Ifremer, AIM and Foselev 

251 Marine for the sea operation and the CC-IN2P3 for the computing facilities. 



252 References 



253 [1] Becker, J.K. 2008, Phys. Rep., 458, 173. 

254 [2] Urry, CM., Padovani, P. 1995, PASP, 107, 803. 

255 [3] Bloom, S.D., Marscher A.P. 1996, ApJ, 461, 657; Maraschi, L., Ghisellini, G., 

256 Celotti, A. 1992, ApJL, 397, L5; Dermer CD., Schlickeiser R. 1993, ApJ, 416, 

257 458; Sikora, M., Begelman, M. C, Rees, M. J. 1994, ApJ, 421, 153. 

258 [4] Gaisser, T.K., Halzen, F., Stanev, T. Phys. Rep. 258 (1995) 173; Learned, J.G., 

259 Mannheim, K. Ann. Rev. Nucl. Part. Sci. 50 (2000) 679; Halzen, F., Hooper, D. 

260 Rep. Prog. Phys. 65 (2002) 1025; Neronov, A., Ribordy, M. 2009, Phys.Rev., 

261 D80, 083008; B6ttcher, M. 2007, astrophys. space Sci., 309, 95. 

262 [5] Abdo, A. A. et al. 2010, ApJ, 722, 520. 

263 [6] Ageron, M. et al, Nucl. Instrum. Meth. A656 (2011) 11. 

264 [7] Amram, P. et al, Nucl. Instrum. Meth. A484 (2002) 369. 

265 [8] Aguilar, J. A. et al, Nucl. Instrum. Meth. A570 (2007) 107. 

266 [9] Aguilar, J.A. et al, Astropart. Phys. 34 (2011) 539. 

267 [10] Agrawal, V., Gaisser, T.K., Lipari, P. and Stanev, T. Phys. Rev. D53 (1996), 

268 1314. 

269 [11] Becherini, Y., Margiotta, A., Sioli, M., Spurio, M. Astropart. Phys. 25 (2006) 

270 1; Carminati, G, Bazzotti, M., Margiotta, A., Spurio, M. Comp. Phys. Comm. 

271 179 (2008) 915. 

272 [12] Aguilar, J.A. et al, Astropart. Phys. 23 (2005) 131. 

273 [13] Heijboer, A. 

274 http://antares.in2p3.fr/Publications/thesis/2004/Aart-Heijboer-phd.pdf, PhD 

275 thesis, Universiteit van Amsterdam, Amsterdam, The Netherlands. 

276 [14] Adrian-Martmez, S. et al, ApJL 2011 accepted, astro-ph.HE/1108.0292vl. 

277 [15] Neyman, J. 1937, Phil. Trans. Royal Soc. London, Series A, 236, 333. 

278 [16] Abdo, A. A. et al 2010, ApJS, 188, 405. 

279 [17] Abdo, A. A. et al 2010, ApJ, 715, 429. 



15 



280 [18] Aharonian, F. et al., A&A 430 865-875 (2005). 

2bi [19] Wagner, S. et al, HEAD 2010 27.06. 

282 [20] Albert, J. et al, 2008a, Science, 320, 1752. 

283 [21] Acciari, V.A. et al, 2008 ApJ 684 73-77. 

284 [22] Monitored source list, |http: / / fermi.gsfc.nasa.gov/ssc/data/access /lat / 



16 



